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We solve the ab initio no-core shell model (NCSM) in the complete 8h£l (iV max = 8) basis for 
A = 7 and A = 8 nuclei with two-nucleon and three-nucleon interactions derived within chiral 
effective field theory (EFT). We find that including the chiral EFT three-nucleon interaction in 
the Hamiltonian improves overall agreement with experimental binding energies, excitation spectra, 
transitions and electromagnetic moments. We predict states that exhibit sensitivity to including 
the chiral EFT three-nucleon interaction but are not yet known experimentally. 

PACS numbers: 21.30.Fe, 21.60.Cs, 27.20.+n 



I. INTRODUCTION 

One of the most challenging problems in nuclear 
physics is to calculate nuclear properties starting 
from the strong interactions that accurately describe 
the nucleon-nucleon, three-nucleon and, possibly, four- 
nucleon systems. There are two major issues to over- 
come. First, the basic interactions among nucleons are 
complicated, they are not uniquely denned and there is 
ample evidence that more than just two-nucleon forces 
are important. Second, as a consequence of the com- 
plex nature of the inter-nucleon interactions, the quan- 
tum many-body problem for these strongly-interacting 
self-bound nuclei is very difficult to solve with good pre- 
cision. 

Interactions among nucleons are governed by QCD. In 
the low-energy regime relevant to nuclear structure, QCD 
is non-perturbative and hard to solve directly to obtain 
these inter-nucleon interactions. New theoretical devel- 
opments, however, allow us to connect QCD with low- 
energy nuclear physics through the promising bridge of 
chiral effective field theory (xEFT) jl|. The xEFT that 
includes pions but omits explicit nucleon excitations, pre- 
dicts, along with the nucleon-nucleon (NN) interaction 
at the leading order, a three-nucleon (NNN) interaction 
starting at the 3rd order (ncxt-to-next-to-leading order or 
N 2 LO) [H-[il, and even a four-nucleon (NNNN) interac- 
tion starting at the 4th order (N 3 LO) Q. The details of 
QCD dynamics are contained in parameters, low-energy 
constants (LECs), not fixed by the symmetry. These 
parameters can be constrained by experiment. A cru- 
cial feature of xEFT is the consistency between the NN, 
NNN and NNNN parts. As a consequence, at N 2 LO 
and N 3 LO, except for two LECs, assigned to two NNN 
diagrams, the potential is fully constrained by the pa- 
rameters defining the NN interaction. 

We have previously performed extensive calculations 
for light nuclei with the xEFT NN and NNN interac- 



tions within the ab initio no-core shell model (NCSM) 5]. 
In particular, we investigated A = 6 and 7 nuclei 6], 
A = 6, 10, 11, 12 and 13 nuclei Q, and A = 14 nuclei 
The major conclusion obtained from these calculations 
was the confirmation of the significance of the NNN in- 
teraction not only for the binding energies but also for 
the description of excitation energies and other observ- 
ables such as magnetic dipole (Ml) and Gamow- Teller 
(GT) transitions. The NNN effects were found to be 
enhanced for the mid-p-shcll nuclei 0, where, for ex- 
ample, the 10 B ground-state spin is in agreement with 
experiment only when the NNN interaction is included 
in the Hamiltonian. One of the dramatic consequences 
of including %EFT NNN interactions has recently been 
found essential to explain the anomalous long lifetime 
(suppressed GT matrix element) of 14 C 0. We previ- 
ously discovered another dramatic consequence of NNN 
interactions in producing a strong enhancement of the 
B(M1) transition from the ground state of 12 C to the 
( J 7r , T) = (1 + , 1) excited state Q, a transition that plays 
a major role in inelastic neutrino scattering. This early 
demonstration of the B(M1) enhancement featured the 
use of two realistic NN interactions, Argonne V8' [Toj 
and CD-Bonn 2000 [ll[ each combined with the Tucson- 
Melbourne "prime" NNN interaction This B(M1) 
enhancement has been confirmed with xEFT NN and 
NNN interactions 0]. 

These calculations were performed by employing 
the Okubo-Lee-Suzuki (OLS) effective interaction ap- 
proach [l3|, [3] primarily in the 6HQ (N max = 6) ba- 
sis space. The exceptions are the A = 6 and A = 14 
results which were obtained in the 8M1 space. It is 
desirable to extend all the calculations to larger ba- 
sis sizes for several reasons. First, one would like 
to check the convergence of the smaller-space calcula- 
tions. Second, the soft similarity-renormalization-group 
(SRG) evolved interactions are now available including 
the NNN terms [H, [l(| ■ Variational calculations with 
these interactions require bases bigger than just 6M7 to 
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fully establish the systematic trends. Also, the trends 
in results from different renormalization schemes, such 
as OLS and SRG, need to be compared with each other 
to better understand their advantages and drawbacks. 
Third, the importance-truncation approach has been suc- 
cessfully implemented for the NCSM calculations (l7tll8j |. 
That approach requires benchmarking against exact cal- 
culations in the same N ma xMl basis space. Fourth, 
the NCSM has been extended by the resonating group 
method (NCSM/RGM) for the description of nuclear re- 
actions 0. The NCSM/RGM approach relies on the 
SRG interactions and requires basis expansion beyond 
6M1 

We report here our recent technical advances that al- 
low us to perform the full-space 8M7 calculations for all 
the p-shell nuclei. These technical advances include both 
the developments of the configuration interaction code 
MFDn [2(| and the development of the codes that calcu- 
late the NNN matrix elements. In this paper, we present 
results for A = 7, 8 nuclei using the OLS method. Cal- 
culations for other p-shell nuclei using both the OLS and 
the SRG methods are underway and will be reported sep- 
arately. 

In Sect. mi we briefly describe the NCSM approach and 
the technical advances that allowed us to extend the basis 
size of the calculations. Results for A = 7 and A = 8 
nuclei are given in Sect. Mil Conclusions are drawn in 
Sect. [TV] 



II. AB INITIO NO-CORE SHELL MODEL 

In the ab initio NCSM, we consider a system of A 
point-like non-relativistic nucleons that interact by re- 
alistic NN or NN + NNN interactions. Unlike in stan- 
dard shell model calculations, in the NCSM there is no 
inert core, all the nucleons are considered active - there- 
fore the "no-core" in the name of the approach. Be- 
sides the employment of realistic NN or NN + NNN 
interactions, two other major features characterize the 
NCSM: i) the use of an harmonic oscillator (HO) ba- 
sis truncated by a chosen maximal HO excitation energy 
iV max /jO (cquivalcntly, the number of HO quanta iV max ) 
above the unperturbed ground state (i.e. the lowest pos- 
sible HO configuration) of the A-nucleon system; and 
ii) the use of effective interactions. The reason behind 
the choice of the HO basis is the fact that this is the 
only basis known (aside from the plane wave basis) that 
allows one to use single-nucleon coordinates and conse- 
quently the second-quantization representation without 
violating the translational invariance of the system. The 
powerful techniques based on second quantization and 
developed for standard shell model calculations can then 
be utilized - therefore the "shell model" in the name of 
the approach. As a downside, one has to face the conse- 
quences of the incorrect asymptotic behavior of the HO 
basis. The preservation of translational invariance is a 
consequence of the N^axMl truncation. 



In order obtain a reasonable approximation in a finite 
basis space (characterized by A max and Ml) to the ex- 
act results in a complete (but infinite-dimensional) basis 
space, we construct an OLS effective interaction from 
the original realistic NN or NN + NNN potentials by 
means of a unitary transformation. We carry out this 
transformation at the two-body level ( "NN only" ) and at 
the level including both NN and NNN interactions ( "NN 
+ NNN"). In principle, one can perform the unitary 
transformation and generate many-body interactions up 
to and including all nucleons. However, going beyond NN 
+ NNN is technically very challenging. We may refer to 
our implementation of OLS as the cluster-truncated OLS 
approach. 

The OLS effective interaction depends on the basis pa- 
rameters (JV max ,f!) and recovers the original realistic TV TV 
or NN + NNN interaction as TV max approaches infinity. 
In principle, one can also perform calculations with the 
unmodified, "bare" , original interactions. Such calcula- 
tions are then variational with respect to JV max and Q. 

In this work, we use NN and NNN interactions de- 
rived within the chiral EFT. In particular, we employ 
the chiral N 3 LO NN interactions from Ref. [H [H| and 
the chiral N 2 LO NNN interaction Q in the local form 
of Ref. [23[ . For the low-energy constants of the NNN 
interaction not fixed by the two-nucleon data, we adopt 
values that reproduce the triton binding energy and half 
life J24[. Next, we calculate three-body effective interac- 
tion from the chiral NN + NNN interactions using the 
OLS procedure. As mentioned above, we adopt a cluster 
truncation which means that the three-body interaction 
is derived from full-space three-nucleon system solutions 
and the resulting three-body effective interaction is then 
input into the shell model code for the A-nucleon sys- 
tem. A large-scale diagonalization is then performed in 
the ^4-nucleon N max Ml HO basis. 

As the three-body effective interactions are derived in 
the Jacobi-coordinatc HO basis but the p-shell nuclei cal- 
culations are performed using the shell model code in 
a Cartesian-coordinate single-particle Slater-determinant 
M-scheme basis, we need to perform a suitable trans- 
formation of the interactions. This transformation is a 
generalization of the well-known transformation on the 
two-body level that depends on HO Brody-Moshinsky 
brackets. Details of this transformations are given in 
Refs. [H, [26J. In this work, we use the particular 
version given in the appendix of Ref. @. The corre- 
sponding computer code was improved compared to ear- 
lier applications Q that allow us now to perform the 
transformations up to TV max = 8 basis spaces for all p- 
shell nuclei. We note that for the p-shell nuclei with 
A > 7, the number of NNN M-scheme matrix elements 
increases compared to the A = 6 case that was handled 
up to iVmax = 8 already in Ref. 0. We note that the 
Jacobi-to-Slater-determinant NNN transformation was 
further improved recently by utilizing a factorization with 
a NNN coupled- JT scheme [27| allowing to reach still 
larger 7V max spaces. 
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FIG. 1: (Color online) Calculated ground state energy of 7 Li 
in the NCSM with chiral EFT NN and NNN interactions 
that reproduce the triton binding energy and half life. The 
dependence on the HO frequency and size of the basis is pre- 
sented. 

It is a challenge to utilize the M-scheme NNN inter- 
action in a shell model code. First, one has to deal with 
a large number of the NNN matrix elements and, sec- 
ond, the number of the non-zero many-body Hamilto- 
nian matrix elements increases by more than an order of 
magnitude compared to calculations with just NN in- 
teractions. Both these issues were successfully addressed 
in the newer versions of the shell model code MFDn, a 
hybrid OpenMP and MPI code (20| . The calculations 
discussed in this paper were performed on Franklin at 
NERSC using up to six thousand cores for the largest 
runs. Current versions of MFDn are capable of handling 
dimensions exceeding 1 billion with NN + NNN Hamil- 
tonians. This is sufficient for the N max = 8 basis space 
of any p-shell nucleus. The largest calculations to date 
are the iV max = 8 results for 14 C and 14 N presented in 
Rcf. 0. 

III. APPLICATIONS TO A = l AND A = 8 
NUCLEI 

Our calculations for both A = 7 and A = 8 nuclei were 
performed in model spaces up to 8MI for a wide range of 
HO frequencies. We then selected the HO frequency that 
corresponds to the ground-state energy minimum in the 
8M) space for detailed comparison of our results with ex- 
perimental data. To further elucidate convergence prop- 
erties, we discuss the dependencies of some observables 
on the basis space parameters (N max ,Cl). 

The 7 Li ground-state energy dependence on the HO 
frequency for different model spaces is shown in Fig. [1] 
for the chiral NN + NNN interaction. We observe that 



FIG. 2: (Color online) Calculated excitation energies of the 
lowest 5 excited states of 7 Li in the NCSM with chiral EFT 
NN and NNN interactions that reproduce the triton binding 
energy and half life. The dependence on the HO frequency is 
presented at the two highest values of basis space truncation 
Af mil i = 6 (dashed lines) and 8 (solid lines). The experimental 
excitation energies are shown on the left for comparison. 

the minimum shifts towards smaller HO frequency as the 
basis space increases. There is also a weak but irregular 
trend towards less dependence on the HO frequency as 
would be expected from reduced roles of the neglected 4- 
body and higher-body interactions. To be specific, Fig. 
[1] shows that the ground state (gs) energy covers a range 
of (2.42, 1.14, 1.33) MeV for N max = (4, 6, 8) in the range 
of Ml depicted. 

The pattern of the iV max = (4, 6) curves in Fig. [TJ is 
very similar to those in Ref. @ where different param- 
eters for the chiral NNN interaction were employed for 
7 Li. There is a shift with the present interaction of about 
1 MeV at -/V max = 6 towards greater binding and towards 
better agreement with experiment when compared with 
the results of Ref. @ with the chiral NNN interaction 
closest (called "3NF-A" in Ref. @ ) to the present case. 

We show in Fig. [2] the dependence of the low-lying 
excited states of 7 Li on the harmonic oscillator energy at 
the two highest values of basis space truncation, N maJC = 
6 and 8. Fig. [5] demonstrates the systematic trend to 
improved independence of Ml with increasing iV max for 
the excitation spectra (slopes of the excitation energies 
decrease with increasing iV max ). Furthermore, the shifts 
in the excitation energies when proceeding from N max = 
6 to 8 are less than the spread in the Ml dependence at 
]V max = 6 over the range of Ml depicted in Fig. [2] 

Quantifying the uncertainties in our results for nuclear 
observables is a major challenge. The systematic uncer- 
tainties due to lack of complete convergence dominate our 
overall uncertainties by at least an order of magnitude. 
The actual uncertainties are dependent on the specific ob- 
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scrvable as well as whether the observable is hindered or 
enhanced compared to phenomenological single-particle 
values. To give an example, we estimate the numeri- 
cal uncertainties in our calculated gs energies and excita- 
tion energies to be around (50, 1) keV arising respectively 
from the numerical evaluation of (1) the effective NN or 
NN+NNN interactions for the selected basis space and 
(2) the numerical solution of the many-body eigenvalue 
problem. However, as we shall discuss now, our overall 
uncertainties are dominated by incomplete convergence 
with increasing basis size. Since we do not have con- 
verged calculations to calibrate these uncertainties, we 
will simply quantify specific measures of our basis space 
dependence. Convergence arises when these measures of 
basis space dependence vanish. 

The results in Fig. [2] provide indications of our basis 
space dependence. For the present work, we adopt the 
following procedure to estimate the dependence of the 
excitation energies on the basis-space truncation. We 
quote two quantities: (1) one half of the total spread in 
the excitation energy over the range in Ml shown in this 
figure at N max = 8, and (2) the total shift in excitation 
energy obtained from the increment of N max — 6 to 8 
at the selected optimum frequency. These quantities are 
quoted in parenthesis beside each excitation energy result 
for the NN + NNN interaction in the tables below. For 
the excited states in the present work we present both 
estimates in the above respective order to show the state- 
by-state fluctuations in each quantity. For uniformity 
and completeness, we present two significant figures for 
each estimate. 

Since many of the states we investigate are resonances 
and our basis lacks explicit coupling to the continuum, 
we expect, in accordance with the findings of Ref. [28| . 
that broader resonances are associated with larger Ml- 
dependence in the HO basis calculations. Thus, one may 
also interpret our first measure of basis-space dependence 
as a rough indicator of the resonance width. This may 
be useful for estimating relative widths [28| . 

From the N max = 8 curve in Fig. [T]we select the opti- 
mal frequency as Ml = 13 MeV for examining our results 
in greater detail. This adoption sets one of the inputs to 
the determination of the basis-space dependence in ex- 
citation energies as just described. We also define the 
basis-space dependence of our total gs energy as the dif- 
ference in total energy at this adopted minimum for the 
basis space increment from iV max = 6 to 8. As an ex- 
ample, this produces the estimate of 0.44 MeV for the 
7 Li gs energy which is quoted in parenthesis next to the 
eigenvalue in Table U 

We observe a similarity in the N max dependence or 
results in Figs. Q]and[2J In both cases, our estimated un- 
certainties range up to several hundred keV (see Table IJ). 
However, in the absence of a firm trend in N max for our 
results, one should not take our quoted uncertainties as 
estimates of numerical accuracy but rather as character- 
istics of the dependence of the results on the presently 
available basis spaces. 
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FIG. 3: (Color online) Calculated and experimental excita- 
tion energies of 7 Li. Dependence on the size of the basis is 
presented. The chiral EFT NN and NNN interaction was 
used. The isospin of the states is T — 1/2 unless shown oth- 
erwise. See the text for further details. 



We show the low-lying spectra of 7 Li in Fig. [3] at the 
optimum frequency and at the sequence of iV max trunca- 
tions corresponding to the curves in Fig. Q] The energies, 
radii and electromagnetic observables are summarized in 
Tables [Q and M where we also include the 7 Bc results. 
We obtain the same level ordering for 7 Be and 7 Li which 
is also the same for both NN and the NN + NNN in- 
teractions with the exception of a reversal of the 7/2^ 
and 3/2^ levels in 7 Be. That is, in 7 Be, the experimen- 
tal 7/2^ and 3/2^~ levels are reversed compared to our 
results and the situation in 7 Li. On the other hand, our 
TV TV + NNN ordering is in agreement with experiment 
for the 9 lowest states in 7 Li. 

Our calculated spectra for both of the A = 7 nuclei 
show a reasonable stability with respect to the frequency 
change. The results in Table U (and ^4 = 8 results in 
Tables HTT1 and IVT1 below) indicate that there are residual 
differences between theoretical and experimental energies 
that are significantly larger than our quoted basis-space 
dependence of the calculated results. It will be inter- 
esting to sec if the differences between theory and ex- 
periment persist once more accurate calculations become 
feasible. If they do, the question becomes whether these 
differences are significantly reduced, for example, when 
a chiral NNN interaction becomes available that is more 
complete than the one currently available [2i|. 

We present in Table |TT] a selection of results for mag- 
netic moments, Ml transitions and other properties of 
the A = 7 nuclei. All electromagnetic observables are 
evaluated with the free-space electromagnetic coupling 
constants. That is, we do not employ effective charges or 
effective magnetic moments for the nucleons. 

The results in TableHJwith NN alone and NN+NNN 
interactions are both in reasonable agreement with ex- 
periment. One observes that there is a trend for radii 
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TABLE I: The 7 Be and 7 Li ground and excited state energies 
(in MeV) obtained using the chiral AW and chiral NN + 
NNN interactions. The HO frequency of Ml = 13 MeV and 
the 8hQ model space were used. Our measures of basis-space 
dependence are given for the last two significant figures of 
the quoted theory result. Two quantities, as explained in the 
text, are quoted in parenthesis for excitation energies with the 
notation: (0.5 x total range of swing with Ml at N max — 8; 
difference at htt — 13 MeV between Nmax = 6 and 8 results). 
Only the second quantity is quoted for the magnitude of the 
total ground state energy. The 7 Be states labelled "mixed iso" 
have large isospin mixing and their basis-space dependence 
can be approximated by the dependences in the corresponding 
states of 7 Li. Experimental values are from Ref. [30l |. 
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and quadrupolc moments to increase with increasing ba- 
sis size and/or decreasing frequency. This is, in part, a 
consequence of the incorrect asymptotics of the HO basis 
and also our basis space truncation. Convergence rates 
for the radii and quadrupole moment appear improved 
at a smaller HO frequency since they are less dependent 
on the basis truncation. 

We have performed various tests to establish that our 
calculated electroweak observables, those near or greater 
than single-particle values, are accurate to three signif- 
icant digits. However, we are not able to quantify the 
basis-space dependence for our electroweak observables 
at the present time. By presenting results at different 
values of Ml and N max we provide a preliminary indica- 
tion of those dependences. Quantifying their systematic 



uncertainties more rigorously will require an extensive 
separate investigation. In the meantime, our best esti- 
mate of the exact theoretical result is the one obtained 
in the largest basis space with the optimum value of htt. 

The magnetic moments and B(M1) values tend to be 
about 10% to 20% smaller in magnitude than experi- 
ment. We believe this is an acceptable range of difference 
since we have not included exchange current corrections 
which can easily be in the range needed to improve the 
agreement with experiment. Future work will address 
these corrections. In the meantime, we can offer sup- 
port for this belief by citing the corrections due to two- 
body currents obtained in the ab initio evaluation of the 
the magnetic moment of the 7 Li gs using GFMC tech 
niques with AV18 plus Illinois-2 three-body forces 31 



In that investigation, the two-body currents raised the 
7 Li gs magnetic moment from 2.9 /in to 3.2 fipf. Simi- 
larly, the two-body currents changed the 7 Be gs magnetic 
moment from -1.06 to -1.49 ^n- Both these correc- 
tions are in the direction and of the magnitude needed 
to explain the difference of our results from experiment. 

We adopt the experimental gs rms charge radius for Li 
(2.44(4) fm) and 8 Li (2.34(5) fm) from a recent detailed 
analysis of the experimental and theoretical gs properties 
of the Lithium isotopes [32[ . We also adopt the correc- 
tions they define and evaluate (finite proton charge den- 
sity, neutron charge density, etc.) in order to extract a 
point proton rms radius r p from the measured rms charge 
radius that we quote in our tables as the experimental 
value for r p to be compared with our theoretical results. 
Note that our theoretical r p is free of spurious cm motion 
effects. 

We present our 8 B and 8 Li ground-state and excited- 
state energy results in Table IIIII The basis size de- 
pendence of the 8 B spectra calculated using the chiral 
NN + NNN interaction and the optimal HO frequency 
of Ml = 13 MeV is shown in Fig. |4j Similar conclusions 
can be drawn as for the A = 7 nuclei concerning con- 
vergence. The dependence of the gs energy on the basis 
size and the HO frequency is somewhat larger than was 
observed for the A = 7 nuclei. This may be due, in part, 
to greater proximity to breakup thresholds in the A = 8 
nuclei we investigate here. 

One noticeable difference between the chiral NN and 
the chiral NN + NNN predictions appears among the 
low-lying levels - an interchange of the order of the Of 
and 3* states. We note that the Of state has not been ob- 
served experimentally. However, the recent Ref. [33| does 
claim observation of the low-lying + resonance based on 
the R-matrix analysis of the p- 7 Be scattering experiment 
performed in the energy range between 1.6 to 2.8 MeV 
in the center of mass. They suggest the Of resonance is 
at 1.9 MeV which places it below the experimental 3^ 
state. However, our calculated 0+ energy obtained with 
the chiral NN and NNN interaction is above our calcu- 
lated 3^~ state. On the other hand, note that this state 
has a larger Ml dependence than lower-lying states which 
suggests a proper scattering treatment is needed for its 
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TABLE II: The 7 Be and 7 Li point-proton rms radii (in fm), 
ground-state quadrupole (in efm 2 ), magnetic moments (in 
/xn) and Ml transitions (in fj,^) obtained within the NCSM 
for different HO frequencies (given in MeV) and model spaces 
for the chiral AW and chiral AA + AAA interactions. Most 
experimental values are from Ref. [3(il |. The point proton 
rms radius r p for 7 Be is evaluated from the experimental rms 
charge radius of 2.647(17) fm from Ref [34T | using corrections 
of Ref [32] as discussed in the text. The point proton rms ra- 
dius r p for 7 Li is evaluated from the experimental rms charge 
radius of Ref [33| as discussed in the text. 
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3.71(48) 




7 Li 






chiral AA 






Ml A m ax 




Q 




B(M1; i" 


2 ") 


13 4 


2.130 


-2.563 


3.038 


4.268 




13 6 


2.140 


-2.786 


3.019 


4.178 




13 8 


2.176 


-2.987 


3.003 


4.100 




Expt. 


2.32(5) 


-4.06(8) 


+3.256 


4.92(25) 








chiral AA + AAA 






Ml Amax 




Q 




B(M1; i" ->• 


?. > 


11 4 


2.225 


-2.683 


3.035 


4.324 




11 6 


2.189 


-2.811 


3.023 


4.230 




11 8 


2.187 


-2.936 


3.009 


4.144 




13 4 


2.091 


-2.422 


3.034 


4.260 




13 6 


2.086 


-2.587 


3.012 


4.154 




13 8 


2.114 


-2.752 


2.993 


4.068 




15 4 


2.002 


-2.252 


3.028 


4.188 




15 6 


2.030 


-2.449 


2.995 


4.088 




Expt. 


2.32(5) 


-4.06(8) 


+3.256 


4.92(25) 





TABLE III: The B and Li ground and excited state energies 
(in MeV) obtained using the chiral AA and chiral AA + 
NNN interactions. The HO frequency of Ml = 13 MeV and 
the 8hQ model space were used. See the caption to table U 
and the text for explanation of the basis-space dependences 
quoted in parenthesis. Experimental values are from Ref. [3f|. 



8 B 

Expt. A A A A + AAA 



|£ gs (2+l)| 


37.7378(11) 


31.38 


36.35 (67) 




0.7695(25) 


0.81 


0.95 (16;04) 


£ x (3+l) 


2.32(20) 


2.83 


2.73 (15;09) 


£ x (0+l) 




2.29 


3.70 (80;25) 


E x (l+1) 




3.11 


4.44 (82;27) 


£x(2+l) 




3.66 


4.62 (44;15) 


£ x (2 3 +l) 




5.11 


5.79 (33;22) 


£ X (L+1) 




4.64 


5.85 (66;25) 


Ex (4+1) 




6.27 


7.20 (34;18) 


£ x (3+l) 




7.12 


7.98 (47;26) 


£ x (0+2) 


10.619(9) 


11.15 


11.68 (27;30) 




8 Li 

Expt. 


NN 


iVA + AAA 


1^(2+1)1 


41.277 


34.86 


39.95 (69) 


£x(lfl) 


0.981 


0.86 


1.00 (16;03) 


£ x (3+l) 


2.255(3) 


2.86 


2.75 (16;09) 


£ x (0+l) 




2.51 


4.01 (84;20) 


£ x (l+l) 


3.210 


3.33 


4.73 (84;21) 


£ x (2+l) 




3.78 


4.78 (44;12) 


£ x (2 3 +l) 




5.22 


5.94 (37;20) 


£ x (l+l) 


5.400 


4.81 


6.09 (70;22) 


£x(4+l) 


6.53(20) 


6.44 


7.45 (36;15) 


S x (3+1) 




7.31 


8.24 (50;22) 


£ x (0+2) 


10.822 


11.25 


11.77 (27;29) 



properties. It is known that the positions of resonances 
are affected by the coupling to the continuum as demon- 
strated, e.g. in Ref. [3(| where the 8 B and 8 Li resonances 
were investigated within the NCSM/RGM approach. 

In Table IIIII we also predict a significant number of 
additional levels in these A = 8 systems that are not 
yet known experimentally. The ordering of these pre- 
dicted levels is sensitive to the presence of the NNN in- 
teraction so it would be very valuable to have additional 
experimental information on these states. We also note 
that many of these predicted states have larger basis- 
space dependences which are dominated by their HO fre- 
quency dependence in the N max = 8 basis space. This 
suggests that these continuum states may be somewhat 
broader resonances than the established states since in- 
creasing frequency dependence in HO basis treatments of 
resonances has been correlated with increasing resonance 
width [28[ as mentioned above. 

In Table HVl we compare A = 8 experimental and the- 
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Exp 
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FIG. 4: (Color online) Calculated and experimental excita- 
tion energies of 8 B. Dependence on the size of the basis is 
presented. The chiral EFT AW and NNN interaction was 
used. The isospin of the states is T = 1 unless shown other- 
wise. See the text for further details. 



oretical results for a selection of electromagnetic observ- 
ables as we did above for the A = 7 nuclei. Here, the 
radii and quadrupole moments are somewhat larger and 
closer to experiment in our chiral AWcalculations due, 
in part, to weaker binding. In addition, contrary to the 
7 Be- 7 Li case, we observe here an interesting difference 
between the NN and NN + NNN cases for the mag- 
netic moment prediction. By including the NNN inter- 
action the magnetic moment of 8 Li is significantly greater 
than that of 8 B in agreement with experiment, while the 
NN interaction alone predicts the opposite. Clearly, our 
results suggest that the A = 8 magnetic moments are 
sensitive to a presence of three-nuclcon interaction in the 
Hamiltonian. These A = 8 magnetic moments also show 
sensitivity to Ml since they change by 10% - 20% in the 
largest basis space over the range of Ml values in Table 

El 

Next, consider the B(M1;1 + — > 2 + ) transitions pre- 
sented in Table IIVI Here, the calculated matrix elements 
are 20% - 35% smaller than the experimental values. 
More significantly, the calculated results are nearly un- 
changed when the NNN interaction is included. Both 

these features are reminiscent of the B(M1;| — > § ) 
transitions in the A = 7 nuclei shown in Table |TIJ 

Our calculated gs energies of A = 8 nuclei are summa- 
rized in Table fVl and shown in Fig. [3] With both the NN 
interaction alone and the NN + NNN interaction, the re- 
sults appear reasonably stable in going from N ma x = 6 
to N max = 8. The role of the chiral EFT NNN inter- 
action is to shift all calculated gs energies significantly 
closer to the experimental results. There appears to be 
a tendency to undcrbind these nuclei as one moves away 
from the minimum in the valley of stability. 

Our calculated excitation levels of 8 Be are compared 
to experiment in Table I VII and Fig. [5] We note a good 



TABLE IV: The S B and 8 Li point-proton rms radii (in fm), 
ground-state quadrupole (in efm 2 ) and magnetic (in /xn) mo- 
ments and Ml transitions (in ^) obtained within the NCSM 
for different HO frequencies (given in MeV) and model spaces 
for the chiral NN and chiral NN + NNN interactions. Most 
experimental values are from Ref. [3^]. The point proton 
rms radius r p for 8 Li is evaluated from the experimental rms 
charge radius [12] as discussed in the text. 



8 B 








NN 








Ml iV max 




Q 


M 


B(M1; 1+ -> 


■r 


") 


13 6 


2.409 


5.243 


1.416 


3.175 






13 8 


2.463 


5.569 


1.438 


3.202 






Expt. 




6.83(21) 


1.0355(3) 


4.71(21) 










AW + NNN 








Ml iVmax 


r P 


Q 


A 4 


B(M1; 1+ -> 




■) 


11 6 


2.472 


5.174 


1.296 


3.215 






11 8 


2.484 


5.294 


1.278 


3.253 






13 6 


2.335 


4.464 


1.146 


3.197 






13 8 


2.374 


4.698 


1.157 


3.215 






15 6 


2.253 


4.004 


1.042 


3.187 






15 8 


2.308 


4.317 


1.079 


3.183 






Expt. 




6.83(21) 


1.0355(3) 


4.71(21) 






8 Li 








NN 








Ml iVmax 




Q 


M 


B(M1; 1+ -> 


2" 


■) 


13 4 


2.119 


2.539 


NA 


4.143 






13 6 


2.117 


2.700 


1.272 


4.133 






13 8 


2.147 


2.826 


1.249 


4.125 






Expt. 


2.23(6) +3.27(6) +1.653560(18) 


5.0(1.6) 










AW + NNN 








Ml V m ax 




Q 


A 4 


B(M1; 1+ -> 


2 H 


') 


11 4 


2.221 


2.687 


NA 


4.153 






11 6 


2.177 


2.766 


1.378 


4.186 






11 8 


2.168 


2.830 


1.392 


4.191 






13 4 


2.080 


2.453 


NA 


4.151 






13 6 


2.064 


2.549 


1.500 


4.145 






13 8 


2.085 


2.648 


1.487 


4.127 






15 4 


1.984 


2.245 


NA 


4.172 






15 6 


2.001 


2.394 


1.581 


4.111 






15 8 


2.041 


2.529 


1.543 


4.072 






Expt. 


2.23(6) +3.27(6) +1.653560(18) 


5.0(1.6) 







agreement of the level ordering compared to experiment 
and also a good stability of the spectrum with respect 
to the change of the model space size. An exception in 
this regard is the calculated first excited + state that is 
an intruder state with large multi-/if2 components. The 
appearance of this state and the corresponding 2 + and 4 + 
excitations were discussed in detail in Ref. [33]. We also 
note that the existence of a 2 + intruder broad resonance 
was confirmed in the R-matrix analysis of the reactions 
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Exp 


AW 


NN + NNN 


8 He 


-31.408 


-24.61 


-29.23 (41) 


8 Li 


-41.277 


-34.86 


-39.95 (69) 


8 Be 


-56.499 


-49.70 


-55.33 (84) 


8 B 


-37.737 


-31.38 


-36.35 (67) 


8 C 


-24.782 


-17.86 


-22.22 (33) 



TABLE V: The NCSM results in the 8M1 basis space for the gs 
energies, in MeV, of 8 He, 8 Li, 8 Be, 8 B and 8 C using the chiral 
AW and the chiral AW + NNN interactions. Experimental 
energies are from Ref [38( |. The basis space dependencies, 
explained in the text, are in parenthesis. 
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FIG. 5: (Color online) Calculated and experimental gs ener- 
gies for A — 8. Calculated results are obtained with the chiral 
EFT AW and NN plus NNN interactions at HQ = ISMeV 
and N m ax = 6 and 8. 



with 8 Be as the composite system [39| . 

We note that our present results for 7 Li, 7 Be, 8 Li and 
8 B compare favorably with those of Ref. 14CI obtained 
with the CD-Bonn 2000 and the INOY [4l| NN po- 
tentials. For example, the trend for the magnetic moment 
in B to decrease when one adds the NNN interaction 
is similar to the decrease obtained as one changes from 
the CD-Bonn to the INOY interaction. The magnitudes 
are also comparable. Earlier work [42j with the chiral 
NN interaction alone [U [22|, |43[ investigated the gs en- 
ergy of 7 Li as a function of the scale set in the Similarity 
Renormalization Group (SRC) approach to the effective 
NN interaction and found the gs energy to range from 
-37.8 MeV to -42.0 MeV which spans our own result of 
-38.60(44) MeV in Table |U The dependence on renor- 
malization scale implies that higher-body effective inter- 
actions are required to obtain a stable result. 

Extensive ab initio calculations of A=7 and A=8 nu- 
clei have been performed with Variational Monte-Carlo 
and Green's Function Monte Carlo [44|,[45[ methods using 
NN plus NNN interactions derived from meson-exchange 
theory. These works provide the most precise agreement 
between theory and experiment for the observables they 



TABLE VI: The 8 Be ground and excited state energies (in 
MeV) obtained using the chiral AW and chiral NN + NNN 
interactions. The HO frequency of fifl = 13 MeV and the 8hQ 
model space were used. Experimental values are from Refs. 
[35l [3^|. States flagged with the superscript ( a ) have signifi- 
cant isospin mixing in both the experimental and theoretical 
results. We quote the isospin assigned by Ref. [3^]. Uncer- 
tainties, explained in the text and in the caption to Tabled] 
are in parenthesis. 





Expt. 


NN 


NN + NNN 


1^(0+1)1 


56.499 


49.70 


55.33 (84) 


S x (2+0) 


3.03 


3.45 


3.63 (10;10) 


£x(4+0) 


11.35 


11.87 


12.53 (30;33) 


£ x (2+0 + l) a 


16.63 


16.30 


17.01 (13;18) 


£ x (2+0 + l) a 


16.92 


16.54 


16.77 (12;23) 


E^l+l)" 


17.64 


17.23 


17.91 (07;21) 


£ x (l+0) a 


18.15 


16.87 


18.25 (27;24) 


£*(3+(l)) 


19.07 


19.20 


19.67 (29;25) 


£x(3+(0)) 


19.24 


18.59 


20.07 (25;21) 


£ x (0+l) 




18.78 


20.84 (72;39) 


£x(4+0) 


19.86 




20.47 (50;33) 
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FIG. 6: (Color online) Calculated and experimental excita- 
tion energies of 8 Be. Dependence on the size of the basis is 
presented. The chiral EFT NN and NNN interaction was 
used. The isospin of the states is T — unless shown other- 
wise. See the text for further details. 



investigate and are closer to convergence, in our view. 
In addition, they have evaluated meson-exchange cur- 
rent corrections to clectroweak processes. Though our 
agreement with experiment is not as good and we do 
not yet incorporate exchange current corrections, we are 
able to expand the suite of observables to compare with 
experiment and to provide the platform for systematic 
improvements with anticipated future developments of 
the chiral interactions and exchange current corrections 
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in the chiral approach. 

More recently, a comprehensive review of the Unitary 
Correlation Operator Method (UCOM) [46[ presents ex- 
tensive results for light nuclei including 7 Li using the Ar- 
gonnc V18 NN interaction |47|. With the UCOM method 
they obtain a variational upper bound on the 7 Li ground 
state of -37.4(6) McV which is also close to our result 
though there are many differences in the calculations. 
Their spectra for 7 Li through the first seven excited states 
are also in the experimental order but spread more than 
experiment as are our own spectra. 

Among other efforts with ab initio no-core methods to 
address some of the same nuclei investigated here with, 
we mention our efforts with an NN interaction derived 
by inverse scattering methods, JISP16 [HI applied to 7 Li 
and 8 Li W9| using the no-core full configuration (NCFC) 
method [501 ]. Results of those investigations also appear 
to be in rough accord with the results presented here. 
For example, the gs energy of 7 Li and 8 Li are -38.253(1) 
MeV and -39.485(16) MeV respectively, both within 500 
keV of the results we report here. 

IV. CONCLUSIONS 

In this work, we used the Okubo-Lee-Suzuki renormal- 
ization of the chiral Hamiltonian specific to each model 
space employed and presented results for A = 7 and 
A = 8 nuclei. Our results demonstrate that the NNN in- 
teraction improves the agreement with experimental data 
not only for binding energies but also for excitation en- 
ergies and other observablcs. Among other features, our 
results suggest that the A = 8 magnetic moments are 
especially sensitive to the presence of three-nucleon in- 
teraction in the Hamiltonian. 

Taking into account our estimates of the basis-space 
dependence of our spectra given in Tables U and IIII1 we 
find that there are residual differences between theory 
and experiment that can now be attributed to the need 
for further improvements to our approach. Those im- 
provements could originate from improved chiral 3-body 
interactions, adding chiral 4-body interactions and/or in- 
cluding effective 4-body interactions. We also recall that 
there is imperfect knowledge of the nonperturbative cou- 
pling constants in the currently-employed chiral NN + 
NNN interactions that could, if exploited, remove differ- 
ences between the current theoretical results and exper- 



iment. 

The present results will be useful for comparing with 
calculations performed with the SRG-evolved chiral in- 
teractions that are currently under way (5l| . In this re- 
gard, the extension of the model space to the 8M7 basis 
is significant as already proven in the case of 6 Li calcu- 
lations [ill. 

Wc have now extended the ab initio no-corc shell model 
calculations with two- and thrce-nuclcon forces in the 
complete 8hil basis to all p-shcll nuclei by technical ad- 
vances of the shell model code MFDn and the codes that 
calculate and transform the three-body interaction ma- 
trix elements. As demonstrated in the 8M7 results for 
the A = 14 nuclei Q we have now the capability to 
calculate any p-shcll nucleus in model spaces up to 8hQ 
with matrix dimensions exceeding 1 billion with Hamilto- 
nians that include NNN interactions. This basis exten- 
sion capability is also significant for a further refinement 
of the importance-truncation approach [l8j and for the 
nuclear reaction applications within the NCSM/RGM 
method 

We also note that further improvements in the 
three-body interaction transformation algorithm and 
a new division of tasks between the shell-model and 
the three-body interaction transformation codes will 
allow one to reach even higher iV max fi,f2 basis spaces. 
This is si gnifi cant since both the importance truncation 
approach [H,[13 and the SU3-NCSM [H provide access 
to much higher N max H£l basis spaces with chiral EFT 
interactions. 
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